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Abstract 

We describe a sounding technique that allows us to improve spatial resolution of Irkutsk Incoherent Scat- 
ter Radar without loosing spectral resolution. The technique is based on transmitting of rectangle pulses 
of different duration in various sounding runs and subtracting correlation matrixes. Theoretically and ex- 
perimentally we have shown, that subtraction of the mean-square parameters of the scattered signal for 
different kinds of the sounding signal one from another allows us to solve the problem within the framework 
of quasi-static ionospheric parameters approximation. 



1 Problem 



The problem of the spatial resolution improvement without loosing spectral resolution is one of the basic 
technical problems of the incoherent scatter (IS) technique (as well as other remote sensing techniques, based 
on weak single backscattering from random media). In such techniques the basic parameters that provide 
information about media properties are spectral power, correlation function and correlation matrix of the 
scattered signal (see [1]). One can show that first two functions can be easily obtained from the third one. 
From scattering theory it is also known that correlation matrix of the received signal is an average of the 
irregularities correlation function over the altitudes. The integration over altitudes is made with the weight 
defined by sounding signal shape and processing technique, and statistical averaging < ... > corresponds to 
the averaging over the sounding runs, see [1] : 
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In the literature the integration weight W(r, r) is referred as 'ambiguity function'. The power profile value 
P(r) takes into account the radiowaves propagation, antenna pattern and cross-section dependence on radar 
range, <t(t,t) - correlation characteristics of the scatterers and a(t) - sounding signal shape. 

So, the practical problem of the spatial resolution improvement (formulated as obtaining information 
about correlation characteristics of the scatterers cr(r, t) with integration over the smallest distances range) 
mathematically leads to the problem of the sounding and processing techniques synthesis that gives the 
ambiguity function W(r, r) with necessary characteristics. 

The problem is currently under investigation (see [2]). But, for now, there are two most useful techniques 
for spatial resolution improvement: with using static amplitude-modulated signals (see [3]) and with phase- 
modulated signals, that change theirs shape from one sounding to another (see [4]). Constructing of these 
techniques is valid in the approximation of quasistationary media. In this approximation one can neglect 
changes of the irregularities correlation characteristics during averaging time and averaging over the sounding 
runs becomes equivalent to the statistical averaging. The use of such techniques allows, from one side, to 
obtain high spatial resolution, defined by discrete of temporal variation of the sounding signal, and from the 
other side to obtain correlation characteristics of the irregularities at large delays (lags) , defined by duration 
of the whole signal. 

The well known difficulty of these techniques is impossibility to obtain the correlation function at all 
the lags with the same high spatial resolution and signal-to-noise ratio. This is caused by the following 
peculiarities of these techniques: correlation function of the irregularities at zero lag can not be obtained 



with high spatial resolution; correlation function can not be obtained at all the real lags within pulse duration, 
but only at discrete number of lags; the number of the 'correct' correlation function lags and spatial resolution 
are usually connected: the higher spatial resolution - the more lags in correlation function you get. So one 
can not obtian a lot of lags, necessary for correct estimation of the correlation function, when one need to 
obtain average spatial resolution (about 1/4 - 1/5 of the whole pulse length). 

We present sounding technique that allows to improve spatial resolution of incoherent scatter technique 
without the problems described above. 



2 Solution 

From (JlJ one can see that independence of irregularities correlation function on sounding signal shape 
allows one to use any linear combination of the correlation matrixes obtained during experiment instead 
of simple summarizing over the sounding runs(i.e. averaging over the sounding runs). In this case the 
ambiguity function will be (with accuracy of the custom multiplier) the same linear combination of the 
ambiguity functions for different sounding runs, summarized over the sounding run number j (the fact, for 
example, is used in summation rules of the alternating code technique, see [4]): 
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The ideal ambiguity function must allow to determine irregularities correlation function without distortion 
effects, caused by power profile P(r). At Irkutsk IS radar power profile is very irregular due to both electron 
density profile and Faraday variations profile (see [5]). One can show that the ideal ambiguity function must 
have the shape: 

< W(r, t) >= A(r)B(r) (4) 

We have found that for obtaining the almost ideal ambiguity function we should use a simple rectangle pulse 
as a sounding signal, with constant amplitude and with two custom durations, the pulse duration should 
change from one sounding to another: dom(a 2 „) = [0,Ti ong ]; dom(a2 n +i) = [0,T s hort]; during processing 
we should subtract the correlation matrix for shorter pulse from the correlation matrix for longer pulse. 
In this case Z?2n = l;£*2n+i = — 1 in ([2]). During the subtraction some parts of the ambiguity functions 
are effectively compensate each other, so we call this technique 'effective subtraction technique' (EST). It 
is important to note that (for accurate positioning the ambiguity function) the correlation matrix must be 
calculated in the form: 
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During the experiment, we should calculate correlation matrix ([5]), but separately for sounding by the 
longer pulse Ri on g,2n(T,T) 7 and by shorter pulse R s hort,2n+i(T, t). The correlation matrix Rest{T,t) that 
provides ambiguity function with the almost ideal shape Q is calculated as combination of Riong,2n{T,T) 
and R s hort,2n+l(T,T): 

< ReSt(T,t) >= ^Rlong.2n(T,T) - ^ R s hort,2n+l {T, t) (6) 

The delay T corresponding to the radar range should be measured from the sounding pulse start. Ambiguity 
functions of the correlation matrix ([5| for longer pulse, shorter pulse and the difference between them ^ 
are shown at Fig.l. From the Fig.l one can see that ambiguity function of correlation matrix ^ is close to 
the ideal one @. One can show that the size of ambiguity function over the range (i.e. spatial resolution) 
is proportional to the difference of the pulses duration. Also, one can see, that the size of the ambiguity 
function over the lags (i.e. lag resolution) is proportional to the half-sum of the pulses durations. 
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Figure 1: Ambiguity function for EST 
lation matrix for shorter pulse sounding(B), and for theirs difference (C). Range and lag are in longer pulse 
duration units. 
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Figure 2: The EST correlation matrix 
23/12/2010. 

3 Experimental results 

The EST was tested at Irkutsk IS radar [5] during the experiment 22-23/12/2010. As a sounding pulse 
we used a simple rectangle signal. Its duration was changed during sounding runs and was 900mks for 
longer pulse and 675mks for shorter pulse. The amplitude of the signal remains constant. Using received 
signal, we calculated correlation matrix ^ separately for longer pulse sounding runs Ri on g ,2n{T , t) and for 
shorter pulse R s hort,2n+i(T,T) sounding runs. Correlation matrix that provides high spatial resolution was 
calculated as theirs combination |6]). 

For direct estimation of the ambiguity function we have used sounding runs with observations of the 
localized space scatterers. It can be shown that correlation matrix ^ in this case becomes close to the 
ambiguity function: 

a(r, r) = S(r -Rq);< R ES t(T, t) >=< P(R Q ) >< W(R Q - ~ r) > (7) 

The correlation matrix for localized scatterer (i.e. ambiguity function shape) obtained experimentally for 
EST is shown at Fig. 2. As one can see, new technique actually improves the spatial resolution, as expected 
theoretically. 

At Fig. 3 we present altitude dynamics of the spectral power for daytime and nighttime data, obtained 
with the technique described and with standard technique based only on longer pulse data. From the data 
presented one can see that spectral power valley near zero frequency for the new technique becomes deeper 
and can be observed even at the night data. It should be noted that such a spectral power shape actually 
theoretically expected for the ionospheric conditions (see [1]), so the new technique not only improves spatial 
resolution in comparison with standard technique, but slightly improves spectral resolution too. 
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Figure 3: Spectral power obtained with use of EST (dashed line), and standard Irkutsk IS radar technique 
(solid line) over the same data, obtained during experiment. (A,C,E,G) daytime data, (B,D,F,H) nighttime 
data. Altitudes: (G,H) - 230km, (E,F) - 260km, (C,D) - 290km, (A,B) - 320km 



4 Conclusion 



In conclusion we summarize some properties of the new effective subtraction technique 
resolution is equal for all the lags including zero lag, and correlation matrix can be calculated at any lag with 
almost equal signal-to-noise ratio; the EST uses only long enough simple pulses without any modulation, 
so in the case of the low signal-to-noise ratio (for example at high altitudes) we can use standard signal 
processing technique for both pulses separately with standard spatial resolution, and for high signal-to-noise 
ratio (at low altitudes) - use EST with high spatial resolution; the spatial resolution of EST is defined by 
pulse duration difference, and maximal allowed lags are defined by half-sum of the pulses durations. 
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